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, Abstract. The parton recombination model has turned out to be a valuable 

^ |i tool to describe hadronization in high energy heavy ion collisions. I review the 

• model and revisit recent progress in our understanding of hadron correlations. I 

' also discuss higher Fock states in the hadrons, possible violations of the elliptic 

flow scaling and recombination effects in more dilute systems. 
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1. Introduction 



Hadronization of partons is one of the unsolved problems in Quantum Chromody- 
namics (QCD). Hadronization is a non-perturbative process and deeply connected 
to confinement which by itself is still not understood on a fundamental level. Fortu- 
O I nately, techniques have been invented that help to deal with hadronization phenomena 

2 ' in elementary particle collisions in which hadrons are produced. They are based on 

universal parameterizations of the unknown non-perturbative quantities. A famous 
example are fragmentation functions D^/j^ ^ (0|u|/i(P))(/i(P)|u|0} which parameter- 
k>( \ ize the probability that a hadron h with large momentum P is created from a parton 

}_j ■ u with momentum p > P in the vacuum T. It has been shown that for hadron 

I production at asymptotically large momentum transfer the cross section factorizes 

into a fragmentation function and a perturbative cross section between partons in a 
well-defined way. A second example are exclusive processes where the wave functions 
ip ~ {0\uud\h{P)) of the (few) lowest Fock states of a hadron appear 0. 

Yet another method to deal with hadronization in a simplified way is the quark 
recombination model [H]. As in exclusive hadron production the matrix elements 
tested are of the same form as the one in exclusive processes ~ {0\uud\h). This is only 
a formal correspondence, but attempts have been made to connect recombination to 
both fragmentation W and to exclusive hadron production 5 . 

Recombination can be directly applied to describe hadron production in rather 
dense parton systems 171 I^H 1101 ITTj . The motivation is twofold. First, there is 
an expectation that in a dense parton phase fragmentation is not the correct picture 
for hadron formation. Secondly, it became indeed clear that perturbative hadron 
production and fragmentation are not sufficient to explain RHIC data for hadrons at 
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transverse moraenta of several GeV/c. This conclusion is mainly based on the observed 
large baryon/meson ratio J21j the nuclear modification factors Raa (or Rcp) close 
to or above unity |12[ 1131 [Tl] and the remarkable scaling law for elliptic flow |15[ [H] . 

2. Hadronization of Bulk Matter 

In central heavy ion collisions a hot and dense fireball of deconfined quarks and gluons 
is created. In the recombination model one postulates the existence of thermalized 
parton degrees of freedom at the phase transition temperature Tc which rccombine or 
coalesce into hadrons. It has been found to be sufficient to consider the lowest Fock 
state in each hadron, the valence quarks, which are given constituent masses around 
300 MeV. 

The spectrum of hadrons can be calculated starting from a convolution of Wigner 
functions 7 . For a meson with valence (anti) quarks a and b we have 

(PNm f (PR fd^qd^r^^^ r P ^rP \^ , , 

^-IlJj^J tI^^- - 2' 2 - ^ + 2 ' 2 + 

Here Wab is the 2-particle Wigner function for partons a, b and <I>m is the Wigner 
function of the meson. The sum runs over all possible parton quantum numbers. 
For simplicity the parton Wigner function is usually approximated by a product 
of single particle phase space distributions Wab — WaWb- Several slightly different 
implementations of this formalism have been discussed in the literature [7| |S1 ■ See 
| 10l ITB] for earlier reviews. 

In order to obtain an estimate when recombination is important as a hadronization 
mechanism, one can compare the yields of fragmentation and recombination starting 
from different parton spectra. Thermal parton spectra w ~ e~^l'^ play a special role. 
Recombination of thermal partons leads to an exponential hadron distribution with 
the same slope since 

where x gives the momentum fraction of parton a. Therefore recombination is more 
effective than fragmentation on any thermalized parton ensemble. On the other hand 
one can show that a power-law parton spectrum favors fragmentation at least for large 
Pt, in accordance with perturbative QCD. 

Eq. Q does preserve momentum, but not energy. Therefore it can only be safely 
applied in a kinematic region where mass effects are small, i.e. for Pt S> M . On 
the other hand fragmentation starts to dominate at very high Pt- The large jet 
quenching at RHIC further suppresses the contribution from fragmentation, so that 
recombination effects can be observed at intermediate Pt- 

A very good description of hadron spectra and hadron ratios measured at RHIC 
can be achieved by combining hadron production from recombination for intermediate 
transverse momentum with a perturbative calculation using fragmentation and energy 
loss in the medium 17|. Fig. ^ shows the Pt spectrum of tt'^, p, and A + A in 
central Au+Au collisions obtained in [Jj. The agreement with data is very good for 
Pt > 2 GeV/c. We note that the hadron spectra exhibit an exponential shape up to 
about 4 GeV/c for mesons and up to about 6 GeV/c for baryons, where recombination 
of thermal quarks dominates. Above, the spectra follow a power-law and production 
is dominated by fragmentation. 
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Figure 1. Spectra of tt", p, Kq and A + A as a function of Pt at midrapidity 
in central Au+Au collisions at ^/s = 200 GeV EJ. Dashed lines are hadrons from 
recombination of the thermal phase, dotted line is pQCD with energy loss, solid 
line is the sum of both contributions. Data are from PHENIX (n^, p) I17II12I and 
STAR {K^^, A + A) ISl- 



Let us now assume the parton phase exhibits elhptic flow V2{pt)- Recombination 
makes a prediction for eUiptic flow of any hadron species after recombination [181 17]: 

V2iPT)=nvliPT/n). (3) 

Here n is the number of valence quarks for the hadron. Note that this scahng law 
is derived using the assumption of inflnitcly narrow wave functions, such that the 
momentum is shared equally between the valence partons. Fig. |21 shows the measured 
elliptic flow V2 for several hadron species in a plot with scaled axes V2/n vs Pt/tl. 
All data points (with exception of the pions) fall on one universal curve. This is a an 
impressive confirmation of the quark scaling rule and the recombination model. The 
pions are shifted to lower Py, because most pions in the detectors, even at intermediate 
Pt, are not from hadronization, but from secondary decays of hadrons. This together 
with allowing a finite width of the wave function can improve the description of the 
pion data |T^ . 

The quark scaling of elliptic flow shows that the relevant degrees of freedom at 
early times in the collision are partons and they prove that these partons behave 
collectively. 
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Figure 2. Elliptic flow V2 for 7r+, p, Kq and A as a function of Pt scaled 
by the number of valence quarks n vs Pt/u. The data follows a universal 
curve, impressively confirming the quark scaling law predicted be recombination. 
Deviations for the pions are discussed in the text. Data are taken from PHENIX 
(7r+, p) 15: and STAR (K^, A) 



3. More on Elliptic Flow 

The 4> meson has long been discussed to provide a good test for the vahdity of the 
recombination model |2U| . (jj mesons are as heavy as protons and the questions is 
whether they follow the pattern of the other much lighter mesons, or whether they 
behave like protons and Lambdas. Data from RHIC now impressively confirm that 
the elliptic flow and nuclear modification factors of the 4> are very similar to those for 
kaons ^21,. This is another success for the recombination model. 

Nevertheless one should ask to which accuracy one expects the scaling law for 
elliptic flow to hold. In particular, are the scaling factors of 2 and 3 indeed excluding 
any higher Fock states in the hadrons? A recent study found that higher Fock states 
in an expansion 



could actually be accommodated j221- It is easy to check that the hadron yields 
from a thermal parton spectrum do not change if additional partons are allowed to 
coalesce. Generally speaking the probability to form a cluster on n particles with fixed 
momentum P from a thermal bath is independent of n. 

The situation changes for elliptic flow. Higher Fock states with n partons come 
with their own scaling factor n which seems to destroy the scaling with the number 
of valence quarks. However, under the assumption that the lowest Fock state is 
still dominating, the numerical effect of the corrections is surprisingly small. Fig. 
121 shows the expected violation of the scaling law using the new asymmetry variable 
A — {B ~ M) /{B + M) where B and M are the scaled elliptic flow of a meson and 
a baryon respectively One curve shows A for the lowest Fock state only and the 
other two correspond to hadrons which have a higher Fock state component with one 
additional parton with probability 30% and 50% respectively. Realistic wave functions 



\p) — aoluud) + ai\uudg) + . . . 



(4) 
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Figure 3. Scaling violation variable A = [B — M)/(B + M) for three different 
probabilities C2 to have an additional parton beyond the valence quark structure. 
Wave functions of finite width have been used which lead to a scaling violation 
even in the case C2 = of pure valence quark recombination 12:^1 . 



with finite width have been used which leads to a scahng violation even for a pure 
valence quark configuration. Generally the violations are smaller than 5%. New data 
from STAR analyzes scaling violations in the data and finds them to have the predicted 
sign and order of magnitude |23| . 

The study did not specify the exact nature of the additional partons. They could 
be quark antiquark pairs or gluons. It was discussed for a long time what the fate 
of gluons in the recombination model is. One can argue that they have to dress the 
quarks just above the critical temperature. After this study it is clear that there is 
some room to accommodate gluons during the recombination process directly. Further 
investigations in this direction are necessary. 

4. Hadron Correlations 

More and more data on dihadron correlations are available from the RHIC experiments 
I^I^Sl- The full picture shown by the data seems to be complicated and rich in detail. 
Before one can even begin to analyze the full data set there is one obvious question. 
The measurements find jet-like correlations at intermediate transverse momenta that 
seem to be coming from jet fragmentation rather than recombination. How can this be 
reconciled with the conclusion from single hadron measurements that recombination 
is the dominant source of hadrons in this kinematic regime? 

Originally recombination was successfully applied to describe hadron spectra and 
elliptic flow starting from assumptions about the parton phase at hadronization. One 
crucial simplification always implemented is a factorization of any n-parton Wigner 
function into a product of independent single parton distributions 

n 

Wi,...,n=l[w, (5) 
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By definition this factorization does not permit any correlations between partons. 
Consequently, no hadron correlations can emerge via recombination. It has to be 
emphasized that the above factorization was chosen for simplicity and it was justified 
because single inclusive hadron spectra could be described very well. 

It has been shown in |21] that modifications of (0) including correlations between 
partons indeed lead to correlations between hadrons after recombination. The quality 
of the description of single hadron spectra does not suffer in the process. The source of 
jet-like correlations in the parton medium is the strong coupling of jets to the medium. 
The energy loss is estimated to be up to 14 GeV/fm for a 10 GeV parton OJ]. This 
implies that most jets apart from those close to the surface are completely stopped, 
dumping their energy and momentum into a cell of about 1 fm"^ in the rest frame 
of the medium. This results in a dramatic local heating, creating a hot spot in the 
fireball. Moreover, the directional information of the jet is preserved. Partons of such 
a hot spot exhibit jet-like correlations. 

In a simple extension of the correlation-free factorization ||SJ) was considered 
(see also [H]). It is assumed that correlations are a small effect and that one can 
restrict them to 2-particle correlations Cy . Then a 4-parton Wigner function can be 
written 

Wi234 ^ WiW2W3W4{l + ^Cij). (6) 

i<j 

The correlation functions between parton i and parton j can be arbitrary, but one 
assumes that they vary slowly with momentum and that they are only non-vanishing 
in a subvolume Vc of the fireball. A Gaussian ansatz ~ cqC"*^'^'"'^^^ /i'^'Po) seems to 
be reasonable to describe correlations in azimuthal angle. The 2-meson yield is given 
by a convolution of the partonic Wigner function VF1234 with the Wigner functions ^a, 
$_B of the mesons with an additional integration over the hadronization hypersurface 
E It is assumed that the correlation strength cq ^ 1 which permits omitting 

quadratic terms like Cg or coV2- 

We can now study the associated yield Yab for a given trigger hadron A in a 
given kinematic window as a function of the relative azimuthal angle between the 
two. One finds 

27TNAYABiA<i>) - Qcoe-^^^'>"^^^^°^"NANB. (7) 

The Ni are single particle yields in the kinematic window of the trigger meson or 
associated meson and cq = cqVc/ (tAt) where tAt is the hadronization volume. 
The factor Q = 4 (for two mesons) indicates an enhancement of the correlations 
in the hadron phase compared to the parton phase. The effect is the same as for the 
amplification of elliptic flow by the number n of valence quarks in the hadron. In 
the case of 2-parton correlations, Q counts the number of possible correlated pairs 
between the ua (anti)quarks of meson A and the ub (anti)quarks of meson B. In the 
weak correlation limit where quadratic terms are suppressed only single correlations 
are counted. Apparently one has 

Q = haub, (8) 

thus Q = 6 for a meson-baryon pair and Q — 9 for a baryon-baryon pair. 

Fig.^shows the associated yield of hadrons integrated in azimuthal angle around 
the near side (Acf) — 0) for the case that the trigger is a baryon (proton or antiproton) 
and a meson (pion or kaon) for different centralities. The kinematic window is 2.5 
GeV/c < Pta <4.0 GeV/c for trigger particles and 1.7 GeV/c < Ptb <2.5 GeV/c 
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Figure 4. The associated yield Y^"^" for baryon triggers (right panel) and 
meson triggers (left panel) as a function of Npart- Squares: fragmentation only; 
diamonds: fragmentation and recombination with cq = 0; circles: the same with 
CO = 0.08 X lOO/Nparf, data: PHENIX gS]. 



for associated particles, and IuaI, < 0.35. Fig. 0]shows the associated yield with 
only fragmentation, and fragmentation and recombination both taken into account 
together with PHENIX data [23 . A good description of the data can be reached 
assuming a constant correlation volume. The parameters used for the fireball are the 
same that lead to a good description of single hadron spectra and elliptic flow 21 ■ 

5. Soft-Hard Recombination 

Some implementations of recombination take into account the possibility that hard 
partons can recombine with soft partons |H1 1221 ■ This soft-hard recombination 
smoothens the transition between the pure fragmentation and pure (thermal) 
recombination regions. It was predicted that typical recombination features like a 
larger baryon/meson ratio would then extend to even higher Pt, beyond 6 GeV/c. 
Indeed such behavior might have been seen in recent data |3()j . 

Soft-hard recombination could be a correction to fragmentation which is even 
important in much more dilute systems. The reason is that fragmentation is an 
extremely ineffective mechanism to create baryons. Soft-hard recombination is a very 
good candidate to explain the very larger Cronin enhancement for baryons in d+An 
collisions PJj- There is no thermal parton phase created in d-l-Au collisions, but the 
existence of a cloud of soft partons with exponential spectrum is sufficient to boost 
baryon production via coalescence. Soft-hard recombination could even be responsible 
for the suppression of hadrons at forward rapidities in d+Au (32| . 

Acknowledgments 

I would like to thank B. Miiller, S. A. Bass and C. Nonaka. This work was supported 
by DOE grant DE-FG02-87ER40328. 



Hadronization of Dense Partonic Matter 



8 



References 

Collins J C and Soper D E 1982 Nud. Phys. B 194 445 
Chernyak V L and Zhitnitsky A R 1977 JETP Lett. 25 510 

Efremov A V and Radyushkin A V 1980 Phys. Lett. B 94 245 
Das K P and Hwa R C 1977 Phys. Lett. B 68 459 

Erratum ibid. 73 504 
Migneron R, Jones L M and Lassila K E 1982 Phys. Lett. B 114 189 
Ochiai T 1986 Prog. Theor. Phys. 75 1184 

Fries R J, Miiller B, Nonaka C and Bass S A 2003 Phys. Rev. Lett. 90 202303 

Fries R J, Miiller B, Nonaka C and Bass S A 2003 J. Phys. G: Nucl Part. Phys. 30 223 
Nonaka C, Fries R J and Bass S A 2004 Phys. Lett. B 583 73 
Fries R J, Miiller B, Nonaka C and Bass S A 2003 Phys. Rev. C 68 044902 
Greco V, Ko C M and Levai P 2003 Phys. Rev. Lett. 90 202302 
Greco V, Ko C M and Levai P 2003 Phys. Rev. C 68 034904 
Hwa R C and Yang C B 2003 Phys. Rev. C 67 034902 Hwa R C and Yang C B 2003 Phys. Rev. 
C 67 064902 

Fries R J 2004, J. Phys. G: Nucl. Part. Phys. 30 853 
Fries R J 2005 J. Phys. G: Nucl. Part. Phys. 31 379 

Adcox K et al (PHENIX Collaboration) 2002 Phys. Rev. Lett. 88 242301 
Adler S S et al (PHENIX Collaboration) 2003 Phys. Rev. Lett. 91 172301 
Adler S S et al (PHENIX Collaboration) 2004 Phys. Rev. C 69 034909. 
Long H (STAR Collaboration) 2004 J. Phys. G: Nucl. Part. Phys. 30 193 
Adams J et al (STAR Collaboration) 2004 Phys. Rev. Lett. 92 052302 
Adler S S et al (PHENIX Collaboration) 2003 Phys. Rev. Lett. 91 182301 
Fries R J (2004) Continuous Advances in QCD ed T Gherghetta (World Scientific 2004) 203 

Fries R J 2005 J. Phys. Conf. Ser. 9 290 
Adler S S et al (PHENIX Collaboration) 2003 Phys. Rev. Lett. 91 072301 
Voloshin S A 2003 Nucl. Phys. A 715 379 

Lin Z W and Ko C M 2002 Phys. Rev. Lett. 89 202302 
Molnar D and Voloshin S A 2003 Phys. Rev. Lett. 91 092301 
Greco V and Ko C M 2004 Phys. Rev. C 70 024901 

Nonaka C, Miiller B, Asakawa M, Bass S A and Fries R J (2004) Phys. Rev. G 69 031902 
Pal D (PHENIX Collaboration) 2005 Preprint hcp-ex/0510020 

Cai X Z (STAR Collaboration) 2005 Preprint nucl-ex/05110M] 
MuUer B, Fries R J and Bass S A 2005 Phys. Lett. B 618 77 
borensen P R (STAR Collaboration) 2005 Preprint nucl-ex/0510052] 
Adler C et al (STAR Collaboration) 2003 Phys. Rev. Lett. 90 082302 
Adler S S et al (PHENIX Collaboration) 2005 Phys. Rev. C 71 051902 

Sickles A (PHENIX collaboration) 2004 J. Phys. G: Nucl. Part. Phys. 30 1291 
Fries R J, Bass S A and Miiller B 2005 Phys. Rev. Lett. 94 122301 
Wang X N 2004 Phys. Rev. C 70 031901. 

Fries R J 2005 J. Phys. Conf. Ser. 27 70 

Bass S A, Fries R J and Miiller B 2005 Preprint nncl-t h /Ob 10084| 
Hwa R C and Yang C B 2004 Phys. Rev. C 70 024904 

Hwa R C and Yang C B 2004 Phys. Rev C 70 024905 
Adams J et al (STAR Collaboration) 2006 Preprint nucl- ex/0606003| 
Hwa R C and Yang C B 2004 Phys. Rev. Lett. 93 082302 
Hwa R C and Yang C B 2004 Phys. Rev. C 70 037901 
[32] Hwa R C, Yang C B and Fries R J 2005 Phys. Rev. C 71 024902 
Fries R J 2005 AIP Gonf. Proc. 792 299 



